TO APPEAR IN ApJ LETT 

Preprint typeset using I^I^X style cmulatcapj v. 5/14/03 



oo 
O 
O 



in 



SS433: OBSERVATION OF THE CIRCUMBINARY DISC AND EXTRACTION OF THE SYSTEM MASS 
Katherine M. Blundell 1 , Michael G. Bowler 1 and Linda Schmidtobreick 2 

to appear in ApJ Lett 

ABSTRACT 

The so-called "stationary" Ha line of SS433 is shown to consist of three components. A broad 
component is identified as emitted in that wind from the accretion disc which grows in speed with 
elevation above the plane of the disc. There are two narrow components, one permanently redshiftcd 
and the other permanently to the blue. These are remarkably steady in wavelength and must be 
emitted from a circumbinary ring, orbiting the centre of mass of the system rather than orbiting 
either the compact object or its companion: perhaps the inner rim of an excretion disc. The orbiting 
speed (approximately 200 km s^ 1 ) of this ring material strongly favours a large mass for the enclosed 
system (around 40 M@), a large mass ratio for SS433, a mass for the compact object plus accretion 
disc of ~ 16 Mq and hence the identity of the compact object as a rather massive stellar black hole. 
Subject headings: Stars: Binaries: Close, Stars: Individual: SS433 
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1. INTRODUCTION 

SS433 is famous for the moving emission lines from 
its relativistic processing jets, but the stationary lines 
in the spectrum have enlightened, intr igued and frus 



trated a generation of astronomers ( Crampto n et al 
19801: ICrampton fc Hutchingj Il981|; IP'Odorico eTal 



Its 



19911 iFabrika & Bvchkoval 1 1991 iCies et al.l I2002T ) 
binary nature was first indicated by periodic varia- 
tions in the redshift of stationary H/J and Hel lines 
(|Crampton et al.lll980l ) but the rcdshifts were found to 
be phased incorrectly, relative to eclipses, for these pe- 
riodic shifts to yield reliably the orbital speed. Rather, 
these variations were attributed to H and Hel lines being 
formed in a wind or in gas streams within the system. 
A blue Hell line may be formed at the base of the jets, 
providi ng a measure of the orbital speed of the compact 
object ([F abrika 19971: IFabrika fc Bvchkoval 1990ft : see also 
IFabrika! (120041 

We obt ained a sequence of nightly spectra oyer so me 30 
days (see lBhmdell. Bowler fc Schmidtobreickl2 007) . Our 
spectra did not extend into the blue but we have stud- 
ied in detail the brilliant stationary Balmer Ha emission 
and the daily sampling made obvious properties which 
had escaped earlier observation. Over the period JD 
2453000 +245.5 to +274.5 only one observation is miss- 
ing (+252.5) and in addition SS433 was quiescent and 
very well behaved during this period. After JD +274.5 
observations were more intermittent and at the same 
time SS433 was becoming restive, culminating in flar- 
ing episodes around JD +294.5. This paper is restricted 
to the period terminating on JD +274.5. 

The stationary Ha lines possess a profile more compli- 
cated than a single gaussian and so were fitted, spectrum 
by spectrum, with a number of gaussian components. 
In almost every case three components were required. 
Where wavelengths or redshifts are quoted in this paper, 
they refer to the centroids of the fitted gaussians. 

2. RESULTS 

1 University of Oxford, Department of Physics, Keble Road, Ox- 
ford, OX1 3RH, U.K. 

2 European Southern Observatory, Vitacura, Alonso de Cor- 
dova, Santiago, Chile 



2.1. The wind from the accretion disc 

Fig. 1 shows the components of the stationary Ha line 
as a function of time. The horizontal scale is wavelength, 
time increases upwards and the height of each tick is pro- 
portional to the standard deviation of the fitted gaussian. 
The nightly observations have created a striking pattern 
in this figure: the central component (whose fitted gaus- 
sians have standard deviations as much as 20 A) varies 
in redshift periodically, between two narrow components 
(having standard deviations of only a few A) which dis- 
play no significant motion. The broad component is very 
symmetric which implies it is formed in a fairly trans- 
parent wind; Fig. 2 establishes unambiguously that this 
is so. This figure displays the variation with time of the 
centre of the wind component (Fig. 2a) and of its width, 
measured by the fitted standard deviation (Fig. 2b) . The 
variation of the central redshift exhibits a periodicity of 
approximately 13 days but the width — a measure of 
the speed of the wind along the line of sight — reveals 
a distinctive pattern with approximately half that peri- 
odicity, which matches the nodding period seen in the 
jets. Fig. 2b also shows that the projected wind speed 
drops as the precessional phase advances (as the jet axes 
lie more in the plane of the sky and the line-of-sight lies 
closer to the plane of the disc). In other data the speed 
of the wind measured primarily in absorption is slow for 
a dense equatorial component but otherwise varies as the 
square of the cosine of the angle between the jets and the 
line of sight (Fabrika 1997; see also Fabrika 2004, pi 14). 
For this reason Fig. 2c is the cosine of the angle between 
the jets of SS433 and the line of sight in th ese same data 
(|Blundell. Bowler fc Schmidtobreickll2007T ) and the ratio 
of the width of the wind line to the square of this co- 
sine is shown in Fig. 2d. In this ratio the 6.5-day pe- 
riodicity has disappeared and the ratio is almost flat, 
apart from a possible 13-day period effect (having max- 
ima at orbi tal phase 0.55, as for the 13-day jet-spe ed 
periodicity ([Blundell. Bowler fc Schmidtobreickl l2007h ). 
Thus the wind increases in speed with elevation above 
the plane of the accretion disc. The dependence of the 
projected speed of the wind (Fig. 2b) on both the pre- 
cessional phase and the nodding of the disc establishes 
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Fig. 2. — Properties of the wind from the accretion disc, (a) 
displays the rcdshift of the wind centroid (grey in Fig. 1) as a func- 
tion of time. A 13-day periodicity is apparent but the redshift docs 
not display a clean cosinusoidal variation. In (b) is displayed the 
standard deviation of the gaussian fitted to the wind; a measure of 
the speed along the line of sight for a transparent wind. There is 
a clear 6.5-day periodicity which matches the nodding of the disc. 
The latter is shown in (c) where cos \ (the cosine of the angle the 
jet axis makes with our line of sight) is plotted, extracted from the 
nodding of the relativistic jets. Finally (d) displays the result of 
dividing the wind speed by the square of the cosine of the angle 
displayed in (c). There may be a residual 13-day oscillation but 
otherwise the wind speed so corrected is flat (but see text), demon- 
strating that the wind is rooted in the accretion disc. The subtle 
grey lines indicate orbital phases of and 0.5; phase occurs at 
Day +254.9. 



Fig. 1. — Upper panel: examples of two spectra of the station- 
ary Balmer H-ct line observed half an orbital period apart. Lower 
panel: wavelengths of components of the gaussians fitted to our 
sequence of spectra. The Julian Date increases vertically and the 
heights of the tick marks are the standard deviations of the fitted 
gaussians. The broad component, in grey, which is identified with 
the wind from the accretion disc, wanders with 13-day period be- 
tween the two almost parallel tracks from the circumbinary disc 
(black). These tracks are separated by almost 400 km s -1 . 



unambiguously that the broad component of the station- 
ary Ha line is formed in the wind from the accretion disc. 

The amplitude of the 13-day oscillation of the centre 
of the wind (Fig. 2a) is < 130 km s _1 and the maximum 
rcdshift occurs just before orbital phase of (eclipse by 
the companion) and certainly after orbital phase 0.75, 



where the compact object is receding most rapidly. It 
would be dangerous to identify this with the speed at 
which the disc is carried round the system, because the 
oscillation is not a good sine curve and the amplitude 
appears to shrink with time. Further, Fig. 2d suggests 
a 13-day periodicity in the wind speed after all the ef- 
fects of orientation of the disc have been taken out. The 
wind from the disc may be interfered with by the larger 
environment. 

2.2. The circumbinary ring 

The redshifts of the two narrow components of Ha are 
remarkably constant. The natural interpretation is that 
we are looking more or less edgc-on to an orbiting ring of 
glowing material and the narrow peaks arc contributed 
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by the two regions tangential to the line of sight. There 
is no significant 13-day period in the sum of the redshifts 
(Fig. 3a): the source must be circumbinary. 

This structure cannot be part of the accretion disc. 
First, the current estimate of the or bital sp e ed of 
the compact object is 175 km s _1 (Fab rikal 119971 : 
iFabrika fe Bvc hkova 1990)- Second, the core of the wind 
oscillates with a redshift amplitude of about 130 km s _1 
and that wind is rooted in the accretion disc (Sec. l2.1|) . 
Finally, in our data after day +274 the accretion disc 
is unveiled and observed to orbit the binary centre of 
mass at approximately 175 km s~ 1 (Blundell, Bowler and 
Schmidtobreick, in preparation). 

The structure from which the pair of fairly stationary 
lines emanate is most unlikely to be a disc orbiting the 
companion donor star. However massive the compan- 
ion may be, it is unlikely to be surrounded by a stable 
excretion disc if it fills its Roche lobe, spilling material 
onto the compact object. Any 13-day periodicity in the 
mean of the disc redshifts, phased as for a disc about 
the companion, is less than 15 km s . A donor orbital 
velocity of 15 km s _1 or less about the binary centre of 
mass would imply that the mass ratio q = m/M (m be- 
ing the mass of the compact object and M that of the 
companion) is less than 0.1 (because q is given by the 
ratio of the orbital speed of the companion to that of the 
compact object). The mass of the compact object would 
then be less than 1 M . 
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Fig. 3. — The upper panel shows the rotational velocity of the 
disc material, obtained from half the difference of the redshifts of 
the red and blue components in Fig. 1 (black). The mean redshift 
is shown as open circles. In the simple case of both tangent regions 
radiating with equal intensity both quantities would be constant 
with time. The lower panel shows the difference between the areas 
of the gaussians fitted to the red and blue components. Blue dom- 
inates near JD +247 and +260; red in the regions of days +255 
and +268. The subtle grey lines indicate orbital phases of and 
0.5; orbital phase occurs at Day +254.9. 



Thus the structure radiating Ha is almost certain to be 
circumbinary and the speed with which this material or- 
bits is easily extracted, for the rather constant redshifts 
and blueshifts shown in Fig. 1. Fig. 3a shows the mean 
redshift and also half the difference, which represents the 



speed at which the material in the glowing ring is rotat- 
ing, under the assumption that the centroids of the fitted 
gaussians correspond to the tangent redshifts. This ap- 
pears to increase slightly with time, from 175 km s -1 at 
precessional phase zero to over 200 km s -1 at precessional 
phase 0.17. This might be due to projection effects as 
the true rotational speed comes closer to the line of sight, 
but another explanation is addition of material at smaller 
radii. We note that the rotational speed in Fig. 3a has ap- 
parent minima at intervals of about 6.5 days. The effect 
is small but not likely to be due to statistical anomalies 
and we remark it is not an effect of the rotating gravita- 
tional quadrupole formed by the compact object and the 
companion. We estimate that this effect would modulate 
the rotational velocity with a period of 6.5 days and an 
amplitude of up to 8 km s -1 , but that the minimum or- 
bital speeds would be observed at orbital phases of 0.25 
and 0.75, rather than appr oximately 0.5 and 1. (See also 
iRudak fe Paczvnskl (|1981l ).) The observed effect is not a 
real change in local disc velocity, but caused by alternat- 
ing intensities in the red and blue components, discussed 
below, and would conceal the tidal effects. 

The systemic velocity is about 70kms _1 but there is 
no indication of whether this is the real recession veloc- 
ity of SS433 or an effect of reddening on the blue edges 
of the lines. A systemic velocity of about 70kms _1 is 
consistent with t he HI velocity o f the environment (see 
recent HI data of lLockman et all (|2007h ). but Hell emis- 
sion data apparent l y prefer a systemic velocit y close to 
zero (jFabrikal 119971 IFabrika fc Bvchkoval fl99"f)h . A n ex- 
otic system such as SS433 may be possessed of a peculiar 
velocity and the true systemic velocity is hard to estab- 
lish through optical spectrosc opy in a filthy environ ment. 
VLBI astrometry reported in lLockman et all (|2007l ) has 
yielded a radial peculiar velocity relative to the local 
standard of rest of SS433 of roughly —17 kms -1 . 

Although these narrow stationary Ha lines stay put, 
their intensities oscillate with a 13-day period. The os- 
cillations of the two are out of phase, such that when 
the redshifted line is most intense, the blue is least 
inte nse. This phenomenon is clear in the raw spec- 
tra ((Schmidto breick fc B lundell 2006) and quantified in 
Fig. 3b, where the difference in areas of the gaussians fit- 
ted to the red and blue lines is plotted as a function of 
time. The inference is of ejected hot gas, cooling and fad- 
ing with time as the source and ejected gas follow their 
respective orbits. This circumbinary stru cture may feed 
the la rge scale feature known as the ruff (Blundell et al. 
l200lh . 

Stationary lines from He I and the Paschen H series 
are an order of magnitude less intense than the bril- 
liant Balmer Ha. They too contain narrow components 
which are radiated by material in the inner circumbinary 
ring and the intensities of the red and blue components 
alternate as the components of Ha. These lines fade 
faster than the Ha components and do not present as 
simple a picture as Fig. 1. We intend to discuss these 
at a later date but here simply mention that both He 
I and the Paschen lines are formed in regions orbiting 
with speeds of approximately 190 km s _1 . The alternat- 
ing intensities of the Paschen lines have been observed 
passing thro ugh the configuration c aptured over a few 
days only by Filippcnko et al. (1988). Both the Hel and 
Paschen lines exhibit much more marked 6.5-day peri- 
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odicities in the difference velocities than Ha, because of 
their stronger alternation in intensities. (Ha and He I 
can b e directly compared in iSchmidtobreick fc Blundell 

3. DISCUSSION AND CONCLUSIONS 

The oscillating redshift of the overall Balmer Ha line 
(|Gies et al.ll2002f) is due to the wind component oscillat- 
ing in redshift and the red and blue narrow components 
alternating in intensity. There may be emissions from 
gas streams within the binary system itself but the dom- 
inant sources of the Balmer Ha emission have now been 
elucidated. 

The glowing ring contributing the narrow lines is orbit- 
ing the system as a whole, rather than either the compact 
object or the donor, and the emission fades with time, the 
hottest spot matching the rotation of the binary system 
and the tail fading, like a katherine wheel. The glowing 
ring is most likely to be located at or close to the inner 
regions of a circumbinary annulus, because the activity 
requires continual refreshment and is linked to the or- 
bital period of the binary. The radiating material might 
either be spilling out piping hot or perhaps disc material 
excited by X-rays from the accreting system. A further 
possibility is that the hot spot migrating round the disc 
is produced wher e excreta hit the i nner r adius of the cir- 
cumbinary ring. iFilippenko et al.l lfl9 88) suggested as a 
possible interpretation of their split Paschen series lines 
an excretion disc fed by overflow of gas from the L2 point 
in the binary system, with the implication the system is 
massive, greater than 40 M Q . Fabrika (1993) remarked 
such an excretion disc should be visible in Ha and the 
intensities of the two lines should vary with the phase 
of the orbit; our observations are in excellent agreement 
with these qualitative predictions. 

The masses of the system and the compact object 
can be extracted from these observations of the cir- 
cumbinary ring and the orbital speed of the com- 
pact object. The radiating material must be orbiting 
close to the inner stable radius for test particles (or 
gas) outside a binary system. This critical distance, 
re-, from the binary centre of mass has been investi- 
gated analytically and also using numerical techniques 
(lArtvmowicz fc Lubowl 119941: IHolman fc Wiegert] 119991 : 
iRudak fc Paczvnskilll981l : ISzebehelvlll980f) . Forbinaries 
with essentially circular orbits the critical radius is 

rc~f(q)a B =.f(q)r(l+q) (1) 

where as is the distance between the components of the 
binary, q is the mass ratio m/M of the components and r 
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is the radius at which the less massive component orbits 
the centre of mass. The value of f(q) is approximately 2; 
more precisely, if f(q) = 2 at q = 0. 1, it has risen to 2.2 
at q = 0.2 and to 2.3 beyond q ~ 0.3 (iHolman fc Wiegert! 
119991 ) . For comparison, the ratio of the critical radius to 
the L2 radius is approximately 1.6, roughly independent 
of the mass ratio q for q > 0.1. 

The speed vc of Keplerian orbits close to the inner 
radius of the circumbinary ring is given by the relation 



v c = {l + q)v x /yJW) (2) 

where vx is the orbital speed of the compact object. 

Thus if f(q) is 2 and vc >vx it follows that the mass 
ratio q > 0.4. The mass of the system is given (for an 
orbital period of 13.08 days) by 

M(l + q) = 1.35/( 9 ) 3/2 («c/100) 3 (3) 
in Mq, vc in km s . 
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Table 1 

Mass of compact object as a function of w and vc- 
Masses are in units of Mq and speeds in kms -1 . 



If the Ha radiation from the circumbinary ring is 
coming from material orbiting at the critical radius or 
greater, then vc is greater than or approximately equal 
to 200 kms -1 . The orbital speed vx is p robably 175 
kms -1 ||Fabrikalll997t iFabrika fc Bvchkovalll990D . Thus 
it seems that q > 0.4 and the mass of the system is 
> 16 Mq and of the compact object > 6M Q , including 
the mass of the accretion disc. These are conservative es- 
timates; Table 1, constructed from equations 2 & 3, gives 
the mass m of the compact object and its disc as a func- 
tion of the orbital speed vx of the compact object and 
the speed of the inner circumbinary disc vc- The smart 
money is on the compact object being a rather massive 
stellar black hole. For vx = 175 km s -1 and Vc= 200 
km s -1 the mass of the compact object is 16 M Q and of 
the companion 22 Mq . 
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